Converging sequences in the ab initio no-core shell model 
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We demonstrate the existence of multiple converging sequences in the ab initio no-core shell model. 
By examining the underlying theory of effective operators, we expose the physical foundations for 
the alternative pathways to convergence. This leads us to propose a revised strategy for evaluating 
effective interactions for A-body calculations in restricted model spaces. We suggest that this 
strategy is particularly useful for applications to nuclear processes in which states of both parities 
are used simultaneously, such as for transition rates. We demonstrate the utility of our strategy 
with large-scale calculations in light nuclei. 
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I. INTRODUCTION 

The ab initio no-core shell model (NCSM) is a method 
to solve the full A-body problem for a system of non- 
relativistic particles that interact by realistic two- plus 
three-body forces. A particular feature of the method is 
the use of effective interactions appropriate for the large, 
but finite, harmonic-oscillator (HO) model spaces em- 
ployed in the calculations. Over the last few years, the 
NCSM method has been established as a very valuable 
tool in nuclear physics, aiming for an exact description 
of nuclear structure startingfrom the fundamental inter- 
action between nucleons [Jill 01 ■ 

The objective of this paper is to highlight the existence 
of multiple converging sequences in the NCSM, and the 
opportunities that this property provides in practical ap- 
plications. In particular, we will discuss and illustrate 
a specific converging sequence that relies on the implicit 
freedom in the choice of model-space cutoff when con- 
structing the cluster-approximated effective interactions 
used in the NCSM. These effective interactions are in- 
tended to take into account effects of configurations out- 
side the model space. In general, the effective interac- 
tion is derived from the underlying realistic inter-nucleon 
potential by a unitary transformation [J, 0, 0, 0, 
The procedure aims to reproduce exactly a subset of 
the eigenvalues to the original Hamiltonian in the finite 
model space. The effective interaction, in principle, be- 
comes an A-body operator but is, in practice, usually ap- 
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proximated at the a-body operator level, where a < A. 
This cluster approximation generates, e.g., a dependence 
on the choice of HO frequency, Ml. However, the con- 
struction guarantees convergence to the exact solution, 
for any value of Ml, as the size of the model space in- 
creases. This guarantee holds independent of the rela- 
tionship between the model space used to evaluate the 
cluster approximation and the model space used to carry 
out the diagonalization of the resulting A-body operator, 
as long as both increase. In addition, we note that for 
any finite model space of dimensionality dp, and in the 
limit a — ► A, we obtain exact solutions for dp states of 
the full problem, with flexibility for the choice of physi- 
cal states subject to certain conditions @. Thus, there 
is an arbitrary number of sequences in the NCSM that 
all converge, in principle, to the same result. 

The implicit freedom in generating approximated ef- 
fective interactions, as mentioned above, will be intro- 
duced in section |TT] in connection to a brief description 
of the NCSM formalism. The strategy that we propose 
to employ for certain applications has several advanta- 
geous properties compared to the traditional strategy for 
evaluating effective interactions. These properties will be 
discussed in section IIII1 while the features and applica- 
bility of the revised strategy will be illustrated in sec- 
tion |lVl The ideas presented in this paper are substanti- 
ated by large-scale calculations for light nuclear systems. 
However, we note that our revised strategy was already 
utilized in two recent papers applying the NCSM to heav- 
ier systems (with A = 47 — 49) in restricted many-body 
model spaces [HI EH ■ Conclusions and perspectives are 
presented in section fVl 
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II. THEORY 

The goal is to solve the A-body Schrodinger equation 
with an intrinsic Hamiltonian Ha = Tr e i+V, where T Te \ is 
the relative kinetic energy and V is the sum of two-body 
nuclear and Coulomb interactions. The NCSM method 
also allows for the inclusion of three-body forces [l2T ]. 
However, while realistic three-nucleon forces have been 
shown to be important in obtaining the nuclear spec- 
tra [13, [H, G3 and in describing electromagnetic and 
weak form factors [15[ , we consider only two-body inter- 
actions in this study. We solve the many-body problem 
in a finite HO basis space and the usual approach is there- 
fore to derive a model-space dependent effective Hamilto- 
nian. For this purpose, we perform a unitary transforma- 
tion of the Hamiltonian, in the spirit of Da Providencia 
and Shakin [|[ and Lee, Suzuki and Okamoto 0, H, 0, [1] , 
which is able to accommodate short-range correlations. 
However, the very first step is to add a center-of-mass 
(CM) HO Hamiltonian to the intrinsic Hamiltonian. In 
the full Hilbcrt space the added H^ M term has no influ- 
ence on the intrinsic properties. When we introduce our 
cluster approximation below, the added Hq M term gen- 
erates a real dependence on the choice of HO frequency, 
but it also facilitates faster convergence to exact results 
with increasing basis size. 

Note that even if the original Hamiltonian contained 
just one- and two-body terms, the transformed Hamilto- 
nian TL contain up to A-body terms. Obtaining the exact 
transformation operator is equivalent to solving the ini- 
tial problem, which would make the procedure imprac- 
tical. Therefore, we introduce the cluster approxima- 
tion. The approximation consists in obtaining the effec- 
tive interaction from the decoupling condition between 
the model space (P) and the excluded space (Q) for the 
a-body problem, where a < A, and then using the ef- 
fective interaction thus obtained in the desired A-body 
problem. See, e.g., Ref. dLULdll for details on the proce- 
dure. This approximation introduces a real dependence 
on the oscillator parameter Ml. The resulting effective 
a-body effective interactions also depends on the nucleoli 
number A and on -/V max , the maximum A-body HO ex- 
citation energy defining the model space. The usual ap- 
proach to handle the ^-dependence is to search for a 
range of Ml values over which the results are weakly /in- 
dependent. This empirical choice then corresponds to the 
optimal HO frequency for the system. 

There are two limiting cases of the cluster approxima- 
tion: First, when a — > A, the solution becomes exact; a 
higher-order cluster is a better approximation and was 
shown to increase the rate of convergence ■ Second, 
when P — > 1, the effective interaction approaches the 
bare interaction; and as a result, the effects of the cluster 
approximation can be minimized by increasing as much 
as possible the size of the model space. 

In this work, we present results obtained at the two- 
body cluster level. In practice, the exact (to numerical 
precision) solutions for the a = 2 cluster are obtained in 



basis spaces of several hundred Ml in each relative mo- 
tion NN channel. The resulting effective Hamiltonian, 
now consisting of a relative two-body operator, and the 
subtracted pure Hq M term introduced earlier, is then in- 
serted into an m-scheme Lanczos diagonalization process 
to obtain the A-body, P-space eigenvalues and eigenvec- 
tors. The evaluation of the A-nucleon Hamiltonian and 
its diagonalization is a highly nontrivial problem due to 
the very large dimensions we encounter. For the present 
work, we performed the many-body calculation with two 
completely independent shell-model codes: a specialized 
version of the code Antoine [13] ; and the many-fermion 
dynamics MFD shell-model code [l8j . At the diagonal- 
ization stage we also add the Lawson projection term, 
(3{Hq M — %Mt) (with (3 being a large positive coefficient) 
to separate the physically interesting states with Os CM 
motion from those with excited CM motion. We retain 
only the eigenstates with pure Os CM motion when eval- 
uating observables. All observables that are expressible 
as functions of relative coordinates, such as the rms ra- 
dius and radial densities, are then evaluated free of CM 
motion effects. 



III. CONVERGING SEQUENCES 

We define our P-space to consist of all A-body configu- 
rations in the oscillator basis with total oscillator energy 
less than or equal some cutoff value (N m + 3A/2)Ml, 
where N m = N m i n + A max , and N m is the sum of 2n + I 
values of the occupied single-particle states in the con- 
figuration. A m ; n is the minimum value required by the 
Pauli principle. As an example, AT m ; n = 2 for 6 Li. The P- 
space is equally described by the cutoff parameter N max , 
that begins with 0. We note that the usual convention 
is to solve only for states whose parity corresponds to 
the configurations in the maximum subspace governed 
by iV max . The opposite-parity states are then obtained 
in the iV max + 1 model space. The corresponding two- 
body cluster model space, P%, is defined by the range of 
two-body states encountered in the P-space. This im- 
plies that the relative (n, I) states are restricted by the 
condition 2n + I < N m - N spsmin , where A spsmin de- 
notes the minimum possible number of HO excitations 
of the (A — 2) spectators. Consider, e.g., a 6 Li calcula- 
tion in the N max = 4 model space. In this case N m = 6 
and iV S p Sm i n = (since the 4 spectator nucleons can all 
fit in the Os shell), which leads to a two-body cutoff at 
2n + l < 6. 

Due to our cluster approximation a dependence of our 
results on A max and on HQ arises. For a fixed cluster size 
a, the smaller the basis space, the larger the dependence 
on Ml. The residual N max and Ml dependences can be 
used to infer the uncertainty in our results arising from 
the neglect of effective many-body interactions. 

The usual NCSM strategy has been to evaluate H s for 
each A-body model space. This strategy leads to the use 
of separate H e g for positive- and negative-parity states. 
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The convergence of, e.g., the energy spectrum of natural- 
parity states is then observed by performing calculations 
in a sequence of model spaces W max = 0, 2,4, . . ., with 
corresponding effective interactions, while unnatural- 
parity states are obtained in the W max = 1,3, 5, . . . se- 
quence. However, the theory guarantees convergence to 
the exact results of any sequence that follows a sim- 
ple rule and brings us to the full ^4-body Hilbert space 
(P — > 1) as the model space is increased. In particu- 
lar, we can consider the following rule for evaluating the 
effective Hamiltonian: 



Amax,cff = W max + iVghift, (1) 



where convergence is obtained for W max — > oo, and where 
IVghift = is the conventional choice. We will instead con- 
sider the combination of choosing Whift = 1 for natural- 
parity states and -/V s hift = for unnatural-parity states. 
This choice implies that we will have the same H e g for 
both positive- and negative-parity states in adjoining 
model spaces, e.g., the JV maX:0 ff = 1 effective Hamiltonian 
will be used in both the iV max = and N max = 1 model 
spaces. The logic for the revised strategy stems from 
four considerations: (1) either strategy will converge to 
the exact result in sufficiently large model spaces; (2) the 
relative position of positive- and negative-parity states 
will converge faster when the same effective Hamiltonian 
is used for both of them; (3) for adjoining spaces in heav- 
ier systems, the predominant sets of pairwise interactions 
are in the same configurations with just one pair at a time 
shifting to the larger space; and (4) for electromagnetic 
transitions between states of opposite parity, the theory 
of the corresponding effective operators will be simpli- 
fied. In addition to these four considerations, the revised 
strategy simplifies our work to compute H a f{ since it is re- 
quired only for every other increment in the basis space, 
such as -/V maXje ff = 1, 3, 5, . . ., to evaluate the converging 
sequence. 

In the next section we will explicitly illustrate the first 
two points by presenting results from large-scale calcula- 
tions for two p-shell nuclear systems: 6 Li and 9 Be. The 
third point was highlighted in two recent papers apply- 
ing the NCSM to systems with A = 47 — 49 in restricted 
many-body model spaces fiol, [Til]. In those papers, the 
authors particularly argued that the bulk of the bind- 
ing should not be altered in proceeding from a Ohfl to 
a 1MI model space in A = 48, suggesting the same H e g 
is preferred. Finally, the fourth consideration is useful 
for applications of the NCSM formalism involving tran- 
sitions between states of different parities; and conse- 
quently for the future description of low-energy reactions 
such as electric dipole radiative capture processes. 



IV. RESULTS 

A. Converging sequences: different JVghift 

We have performed calculations for 6 Li up through 
the 16M7 (N ma .x = 16) model space, and for 9 Be up 
through the 10hQ (iV max = 10) model space. The max- 
imum dimensions of the encountered model spaces are 
d P = 7.9 x 10 8 and d P = 5.4 x 10 8 for 6 Li and 9 Be, 
respectively. 

Most of the results presented in this paper are ob- 
tained using the non-local CD-Bonn 2000 (CDB2k) po- 
tential [l9[, which is a charge-dependent AW interac- 
tion based on one-boson exchange. The off-shell behav- 
ior of the CDB2k interaction differs from local poten- 
tials which leads to larger binding energies in nuclear 
few-body systems. Still, the CDB2k gives undcrbinding 
of nuclear many-body systems as is typical for standard 
high-precision AW interactions. In general, the CDB2k 
potential applied in the NCSM gives a good convergence 
rate and a weak HO frequency dependence. 

A more recent realistic, nonlocal AW interaction, ob- 
tained through an inverse scattering analysis of the NN 
data, is able to provide reasonable binding energies of 
p-shell nuclei and is called the "JISP16" interaction [20| . 
We will use this interaction in Sec. HVBI to illustrate the 
value of exponential fits to converging sequences in the 
NCSM. 

The optimal HO frequency for each isotope was found 
by performing a series of calculations for a number of dif- 
ferent frequencies and model spaces. We searched for the 
region in which the dependence on Ml is minimal; and 
we selected this frequency (from the calculation in the 
largest model space) to use in the more detailed investi- 
gation. In this way we found that the optimal frequencies 
for the CDB2k interaction are Ml = 10 MeV for 6 Li and 
Ml = 12 MeV for 9 Bc (see also Ref. (U). Since we will 
also use the N max + 1 effective Hamiltonian for natural- 
parity states, diagonalized in the A^ max model space, we 
verified by explicit calculations that the same choices of 
optimal frequencies apply. 

Note that the choice of optimal frequency is depen- 
dent on the model space used in the calculations. Con- 
sequently, a larger HO frequency (Ml = 13 MeV) was 
employed in an earlier NCSM study of 6 Li, in which the 
calculations were limited to N max < 10 [|[. 

The ground-state binding energy of 6 Li, as a function 
of increasing model space N max , is presented in Fig. [1] 
The results obtained using the traditional cutoff for eval- 
uating the effective interaction, N max ,eE = W max , is rep- 
resented by black circles, while the results obtained with 
AT ma x,eff = W max + 1 is represented by red squares. A 
number of observations can be made: First, these cal- 
culations demonstrate that the NCSM with the cluster- 
approximated effective interaction is not a variational ap- 
proach. Convergence is not necessarily from above. Sec- 
ond, the difference between the two strategies for com- 
puting the effective interactions is the largest for small 
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FIG. 1: (Color online) Ground-state binding energy of 6 Li 
calculated with the CDB2k Hamiltonian with hQ — 10 MeV. 
Results obtained using two-body effective interactions gener- 
ated with iVshift = and iVshift = 1 are compared. 



model spaces. The choice N max c s = N max + 1 reduces 
the excursions in the renormalized results in the lowest 
basis spaces as a function of iV max . Third, both strategies 
converge to the same value as A^ max — > oo. The binding- 
energy difference between the two calculations is plotted 
on a semi- logarithmic scale in the inset of Fig. [TJ We find 
that the convergence rate towards zero difference is ex- 
ponential. The observed difference in our largest model 
space (iV max = 16) is less than 70 keV (corresponding to 
~ 0.2% of the total binding energy). 

We present the results of a similar comparison for 9 Be 
in Fig. [2J The spectroscopy and various properties of 
this nucleus were extensively studied within the NCSM 
approach in Ref. 2l|. The same observations as for 6 Li 
concerning the convergence rate apply for this case, al- 
though with A = 9 we are slightly further from con- 



verged results. Still, we find an exponential convergence 
rate towards zero difference between the two strategies, 
as shown in the lower panel of Fig. [51 We studied three 
different, low- lying natural-parity states and found the 
same convergence pattern. For the ground state calcu- 
lated in the N max = 10 model space, the difference in 
computed binding energy is 350 keV (< 0.7% of the total 
binding energy). 

As an additional remark, we note that the eigenener- 
gies obtained with the iV max + 1 effective interaction are 
always above the ones obtained with the iV max effective 
interaction. This observation persists for other systems 
that we have studied. An indication of the underlying 
reason for this observation can be obtained from the fol- 
lowing argument: The iV max + 1 effective interaction can 
be diagonalized in the iV max + 1 model space, but we 
choose to diagonalize it in the N max model space. Thus, 
we have a truncation of the "available" basis space so the 
result will be above the one with the same Hamiltonian in 
the full "available" space. This argument does not prove 
though, that it will always be above the results in the 
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FIG. 2: (Color online) Upper panel: Ground-state binding 
energy of 9 Be calculated with the CDB2k Hamiltonian with 
hQ. — 12 MeV. Results obtained using two-body effective 
interactions generated with N^m = and N b mh — 1 are 
compared. Lower panel: Difference between binding energies 
computed with N B hnt — and iVshift = 1 for the first three 
natural-parity states of 9 Be. 



N max basis space with the -/V max effective Hamiltonian. 

Up until now we have demonstrated that the two 
strategies converge to the same value, and that the 
Af max + 1 effective interaction gives better renormalized 
results in the smallest (0M1) model space. We also 
want to show that the relative position of positive- and 
negative-parity states converges faster when they are ob- 
tained by diagonalizing the same effective Hamiltonian. 
This claim is supported by the results presented in Fig. [3] 
In this figure, the combined spectrum of natural- and 
unnatural-parity states of 9 Be is plotted. The two panels 
show the evolution of the spectra with increasing model 
spaces for the two strategies: (left) different effective 
Hamiltonians for each separate model space; and (right) 
the same effective Hamiltonian for each pair of model 
spaces, N max = (0-1), (2 -3}, (4-5), . . .. The eval- 
uated experimental spectrum [22| is shown in the left- 
most column. Indeed, we find that the relative position 
of positive- and negative-parity states is better described 
already for small model spaces with the -/V max +1 effective 
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FIG. 3: (Color online) Spectrum of natural- and unnatural-parity states of 9 Be calculated with the CDB2k Hamiltonian with 
ftfi = 12 MeV. The natural-parity states are obtained using two-body effective interactions generated with iVshift = (left 
panel) or iV s hift = 1 (right panel). 



interaction being used for the natural-parity states. 

Exponential fits to the calculated eigencnergies of the 
first natural- and unnatural-parity states are shown in 
Fig. 2J The two sequences of calculations of £7(3/2^,) 
(with iV s hift = and iVshift = lj respectively) constitute 
an example of a series of converging sequences that con- 
verge to the same exact result. Therefore, we carry out 
a constrained fit of the two data sets. We use exponen- 
tial functions and the constraint for the lit is to have a 
common asymptote to both sequences. The fit produces 
a common asymptote of -51.40 MeV with uncertainty of 
about 200 keV. In carrying out this constrained fit, the 
chi-square weights were based on the local derivative with 
respect to iV max so that results with minimal sensitivity 
to iV max receive greater weight. The iV max = — 2 results 
were excluded from the fit. We attempted alternative 
functional forms, such as a + /?iV~2xj tmt found none to 
be as successful as the constrained exponential fits shown 
in Fig.0] Finally, we present in Fig. 0a plot of the excita- 
tion energy of the first unnatural-parity state in 9 Be, rel- 
ative the natural-parity ground state, as a function of the 
model space size for the two different strategies. These 
results illustrate rather clearly that the relative position 
of natural- and positive-parity states is better described 
in the revised strategy where the same effective Hamilto- 
nian is being used for diagonalization in both iV max and 
iV max + 1 model spaces. The dashed lines correspond to 
the differences between the eigenenergy fits of Fig. 01 The 
common asymptote is E x {\/2^) = 2.10 MeV. The differ- 
ence in excitation energy between the two calculations is 
plotted on a semi-logarithmic scale in the inset of Fig. [5] 



B. Converging sequences: different Ml 

Additional examples of series of converging sequences 
in the ab initio NCSM can be found. In the following we 
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FIG. 4: (Color online) Basis size dependence of the calculated 
£ ( 3/2~ ) and £ ( l/2+) binding energies. The CDB2k Hamil- 
tonian with hi} = 12 MeV is used. The binding energy of 
the (natural-parity) ground state is obtained using two-body 
effective interactions generated with iV s hift = (black circles) 
and iVshift — 1 ( re d squares). See text for details. 



present a series of calculations of the ground state eigen- 
value of 6 Li using the JISP16 realistic NN interaction 
[20j | . Fig. [6] displays the results for the bare interaction 
as a function of iV max up through iV max = 12. In this 
case the series consists of the eight sequences of calcu- 
lations for increasing model spaces performed at differ- 
ent HO frequencies, ranging from M7 = 12.5 — 30 MeV. 
The N mgx = results are suppressed. These calculations 
should all converge to the same result as the dependence 
on the choice of HO frequency disappears with increas- 
ing model space. In addition, the procedure is variational 
due to the use of bare interactions. We carry out a con- 
strained fit using exponential functions for each sequence 
at fixed Ml using only the iV max = 6 through iV max = 12 
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FIG. 5: (Color online) Basis size dependence of the calculated 
E x (l/2~l) excitation energy relative to the lowest negative- 
parity state. The CDB2k Hamiltonian with W, = 12 MeV is 
used. The binding energy of the (natural-parity) ground state 
is obtained using two-body effective interactions generated 
with iVshift = (black circles) and JV s hift = 1 (red squares). 
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FIG. 6: (Color online) Ground-state binding energy of 6 Li 
calculated with the JISP16 Hamiltonian for different HO fre- 
quencies (hQ. = 12.5 — 30 MeV). The curves are extrapolated 
using a constrained exponential fit as described in the text. 



results. The constraint for the fit is to have a common 
asymptote to all sequences. The uniformity of the con- 
vergence is striking and produces a common asymptote 
of -31.33 MeV with uncertainty of about 100 keV. Again, 
the chi-square weights used in the fit were based on the 
local derivative with respect to A max so that results with 
minimal sensitivity to A max receive greater weight. We 
attempted alternative functional forms but found none to 
be as successful as the constrained exponential fits shown 
in Fig. O 



V. CONCLUSION 

The ab initio NCSM is generally characterized by pro- 
viding very fast convergence of many observables with 
increasing model space. This property is obviously very 
valuable when applying the method to studies of many- 
body nuclear systems. The existence of multiple con- 
verging sequences is another important property of the 
method; but one that has not been extensively uti- 
lized in the past. In this paper we have demonstrated 
the existence of multiple converging sequences in the ab 
initio NCSM, and we have discussed some benefits of 
this property for certain applications. In particular, we 
have proposed a revised strategy for computing cluster- 
approximated effective interactions for A-body calcula- 
tions in restricted model spaces. 

The fundamental principle that motivated this work is 
a particular property of the theory of effective operators 
employed in the NCSM; namely that convergence to the 
exact results is guaranteed for any sequence that returns 
the bare Hamiltonian as the model space is increased, 
]V max — - > oo. It was shown in Sees Mil and llVl that the spe- 
cific choice of employing the same effective Hamiltonian 
for calculations in adjacent A-body model spaces resulted 
in some attractive convergence properties. Firstly, it was 
shown that it moderates the excursions in the renormal- 
ized results that are usually encountered when perform- 
ing calculations in smaller model spaces, and secondly, 
that it gives a faster convergence of the relative position 
of binding energies of opposite parity states. We also 
demonstrated, by explicit large-scale calculations, that 
the revised strategy indeed converges to the same result 
as the traditional choice of using a iVshift = effective 
interaction for each model space. We suggested that this 
new strategy is particularly useful for studies of nuclear 
systems in which states of both parities are simultane- 
ously involved. Use of the revised strategy is specifically 
envisioned for calculations of electromagnetic transitions 
between states of opposite parities, as the theory of the 
corresponding effective operator is simplified. 

Another class of converging sequences was demon- 
strated by performing calculations using the bare JISP16 
realistic AW interaction for a series of HO frequencies. In 
this case, our approach obeys the variational principle; 
for each choice of Ml we are guaranteed to approach the 
exact result from above. In addition, all sequences are 
guaranteed to converge to the same result as A max — > oo. 
We utilized this property by carrying out fits to all se- 
quences using exponential functions, with the constraint 
to have a common asymptote. This approach is poten- 
tially very useful for gaining confidence in extrapolated 
results and for estimating their uncertainties. 
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